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Abstract 

Mathematical proofs are presented concerning the existence of solutions of the Maxwell 
equations with suitable boundary conditions. In particular it is stated that the well known 
"delayed potentials" provide effective solutions of the equations, under reasonable conditions 
on the sources of the fields. 

1 Introduction 

Any advanced test on theoretical classical electromagnetism (see for example [1], [2], [3], [4]) 
states that all the physical properties of the electromagnetic fields can be mathematically 
deduced from the Maxwell Equations, and performs such deduction obtaining formulae (the 
"delayed potentials") that allow, in principle, to evaluate the fields, starting from the knowledge 
of the sources, i.e. charges and currents. 

The problem is that the delayed potentials are obtained, as pointed out in the following 
section, with a method which ensures that, if solutions exist for the Maxwell Equations with 
suitable conditions, they must be of the provided form. Therefore the usual deduction is able 
to provide uniqueness theorems, but not existence theorems for the solutions of the Maxwell 
Equations. 

This is not so surprising: also in similar areas, (e. g. the Laplace Equation for harmonic 
functions) the deduction of uniqueness theorems is well easier than the deduction of existence 
theorems. 

For the above reasons, the traditional deduction of the electromagnetism leaves open the 
problem that solutions of the base equations could not exist. The usual way to overcome this 
situation is to state that "due to the physical nature of the problem, solutions must exist" , but 
clearly this cannot be accepted from a mathematical point of view. 

The present paper presents existence theorems, i.e. states that, under suitable conditions on 
the sources, the Maxwell Equations admit solutions, and moreover that the delayed potentials 
provide an effective solution. Similar results could be obtained using the general theory of 
partial differential equations or distribution theory, but the presented deduction is performed 
without any use of such general theories, and requires no mathematical background, besides 
the one required for reading any theoretical test on electromagnetism. 

The first section summarises the standard deduction of the delayed potentials, and points outs 
its defects. 

The second section reports the obtained results, proving that the delayed potentials provide an 
effective solution. 

The third section summarises the results and discusses some extensions. 



2 Standard deduction of electromagnetism 



The Maxwell Equations for the wide space have the form: 
rotE+i^=0 (1) 
'•otH-if = fj (2) 
divH = (3) 
divE = Airp (4) 

where: 



• E is the electric field 

• H is the magnetic field 

• c is the light speed in vacuum 

• j is the current density 

• p is the charge density 

The last two quantities, named "sources of the fields"), are not independent, since they 
must satisfy the continuity equation 

divj + t=0 _ . 

whose physical meaning is the conservation of the total charge. The main problem of the 
electrodynamics is the following: given j and p, as functions of space and time satisfying (5) , 
find the electric and magnetic fields. The standard approach to this problem is summarised in 
the following subsections. 



2.1 Definition of the potentials 

It is trivial to prove that, given any vector function A, the vector function H = rot A satisfies 
(3). Starting from this point, the standard developments of electrodynamics perform the claim 
that, in order to satisfy (3), the magnetic field H must be given by 
H = rotA (6) 

where A is a suitable vector function. A proof of this fact is not difficult, but it is not 
reported, since this result is not required in the following. Assuming that H is given by (6) 
and substituting into (1), the following is obtained: 

rot(E + if^) = (7) 

whose consequence is that the quantity in parenthesis is the gradient of some scalar function 
—ip. Therefore the following formula is obtained: 
E = -gradM-i^ 

Having introduced the scalar potential ip and the vector potential A, the research of the 
fields is restricted to those given by (6) and (7). Consequently (1) and (3) are automatically 
satisfied, while by substituting into (2) and (4) and performing some calculations the following 
two equations are obtained: 

- - grad(divA + 1^) = -^j 

^^-y^ + -cmid^A + l'4) = -4np 
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Normally the last equations are simplified by assuming the Gauge condition: 
divA + it = (8) 

and this assumption is heuristically justified by observing that and A are not uniquely 
defined by (6) and (7). Summarising, the restriction is done to search fields given by (6) and 
(7), where A and (p satisfy (10). With the above assumptions, the equations satisfied by the 
potentials reduce to the form 

AA-^0 = -^j (9) 

A^-^|^ = -47rp (10) 

If solutions of the above equations, subjected to (8), can be found, the conclusion can be 
stated that the Maxwell equations have at least one sohition, given by (6) and (7). 

It is readily seen that (9) and (10) are a system of 4 independent scalar equations; moreover 
they are of the similar form 

AF-^0 = -G (11) 

2.2 Elimination of the time dependence 

Equation (11) is simplified by taking its temporal Fourier transform (accepting that the involved 
functions can be transformed), and by assuming the possibility to exchange integrations with 
differentiations. By standard calculations the following equation is obtained: 
A/ + A:2/ = -ff (12) 

where 



the variable replacing t in the transformed domain is named u (it has a frequency mean- 
ing) 



k = uj/c ^ where c is the light speed in vacuum 
f{oj) and g{ijj) are the Fourier transforms of F and G, respectively. 

In (12) k has simply the function of a constant parameter. 
2.3 Green Theorem 

Equation (12) is treated with the well known Green method, as is done for the Laplace Equation 
for the harmonic functions. As in this case, it is clear that the Green method provides an 
uniqueness theorem for the solution of (14) with assigned values on a boundary, but does not 
provide any existence theorem for such solution. The Green functions is chosen as ^—^ — , where 
d = \/ix — xqY + {y ~ yo)^ + (-^ ~ -s^o)^) and xq, yo, zq are the coordinates of an arbitrary point 
Pq (it is easily verified that this function is everywhere regular and satisfies (12) except at Pq, 
therefore it is an appropriate Green function). 

In such a way the Green Theorem is obtained in the form 

m) = -h I '-^dv + ^ / (/^ - ^-^^/jds 

D S 

and provides the solution of (12) at an arbitrary point interior to a spatial domain D with 
boundary S, using the known term g and the value of the solution and of its normal derivative 
taken on the boundary, in total analogy with the Green formula for the Laplace equation. 
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A similar formula for the solution of (11) is obtained by applying the inverse Fourier trans- 
form, and assuming interchangeability between differentiation and integration: 

D S 

Finally, a behavior at infinity for F is assumed that ensures the vanishing of the second 
integral when the domain G tends to infinity, and the following is obtained: 

F{Po) = -^J^dV (13) 

where presently the integral is extended to the whole space. (13) provides an explicit solution 
of (11), if the assumption is made that a solution exists satisfying all the outlined conditions. 
Coming back to the equations for the potentials, (13) provides: 

A(Po) = -i / i^jM^^zildV (14) 
V>{Po) = -l '^"''''T~^-^ dV (15) 

The above formulae are known as "delayed potentials", and are believed to provide, with 
(6) and (7), a general solution to the electrodynamics problem. 

2.4 Criticism 

As pointed out by the above review of the standard deduction, the fields obtained applying (6) 
and (7) to (14) and(15) provide a solution to the Maxwell equation assuming that the following 
conditions are verified: 

• Existence of a solution which can be derived by the potentials according to (6) and (7) 
(this would not be difficult to prove, if the assumption is made of existence of a solution 
of (1),(2), (3), (4)) 

• Existence of a solution with suitable conditions at infinity 

• As a minor point, validity of various hypothesis concerning the possibility of exchanging 
integrals (singular and with infinite domain) and derivatives 

• Existence of the integrals appearing in (14) and (15) 

• Validity, for the potentials given by (14) and (15), of the condition (8), with is necessary 
to conclude that (9) and (10) provide solutions for the potentials. 

In conclusion, the summarised deduction provides only an heuristic feeling to have found 
possible solutions of the Maxwell equations. 

3 Solubility Theorems 

In order to overcome the summarised limits, the chosen strategy does not attempt to verify the 
validity of the various conditions pointed out at the end of the preceding section, but attacks 
directly the problem to prove that, under suitable conditions on the field sources p and j, (14) 
and (15) provide functions that satisfy (8), (9) and (10), from which it is straightforward to 
prove that the fields obtained by (6) and (7) satisfy the Maxwell equations (1), (2), (3), (4). 
The assumed conditions on p and j are the following: 
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• Regularity, i.e. existence and continuity up to the first derivatives 

• Concerning the space variables, vanishing outside of a finite regular domain D, where 
"regular "means that its boundary admits everywhere a tangent plane 

• Concerning the time dependence, sinusoidal behaviour, or existence of the Fourier trans- 
form, vanishing outside of a finite interval for the transformed variable uj (in practice, 
finitely extended frequency spectrum) 

These conditions, in particular the second and the third, are surely redundant, even if 
acceptable from a physical point of view. The deduction will show some possibility of extensions. 
The proofs are reported for the following cases, of increasing complexity: 

• Electrostatic case 

• Magnetostatic case 

• Monochromatic case 

• General case. 



3.1 Electrostatic case 

This is the case in which all the charges are fixed, i.e. the density p does not depend on the 
time, and the current density j vanishes (note that condition (5) is trivially satisfied). In this 
case the proposed solutions (16) and (17) assume the form 

A(Po) = (14) 
^{P^) = -!l^dV (15) 

D 

where D is the finite domain in which p does not vanish, while the equations to be verified 

((8), (9) and (10)) reduce to 

A<y3 = -47rp (16) 

If a point Pq is external to D , the integral appearing in (15) is not singular and can be 

differentiated under the sign, providing trivially A99 = — 47r/?, as required. If the point Pq is 

internal to D, a transformation to polar coordinates u, a, (3 centered on Pq is performed, giving 
^(Po) = - jT dp J% sin2 a piPo + ^)udu (17) 

in which the vector v = (n sin a cos [3, u sin a sin /3, u cos 0) has been introduced, and g{PQ ,a,f3) 
is the distance between Pq and the point on the boundary of D with polar angles a, j3 , which 
is a regular function. 

In (17) every singularity has disappeared, and this ensures the existence and regularity of 
<p{Pq), and the existence of the vector field given by 

E = — grad((/?) 

The problem is if E can be calculated by taking the derivatives respect to Pq under the 
sign of integral appearing in (19). The response is certainly "yes" if the integrals obtained by 
formal derivation under the sign are uniformly convergent. Performing the formal derivation, 
one obtains 

D 
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where the interrogation mark means that the result must be justified. Performing the same 
transformation to polar coordinates just used, the last integral transforms to 
^(^o) =' - /o"^ d(3 J% sin^ a /,^(^°'"'^) p(Po + v)du 

which is without singularities. Therefore it represents a continuous function of Pq, and the 
uniform convergence is ensured, justifying the formal derivation. It is important to note that 
a similar reasoning could not be applied to the second derivatives, vanishing any attempt to 
verify (20) with a direct calculation. Having justified first order derivatives under the sign of 
ip{Po), the following result is obtained 

E{Po) = Jp{P)j^dV (18) 

with uniform convergence of the integral at its unique singular point Pq. Now let V be an 
arbitrary regular domain internal to D, with boundary S. The aim is to calculate the flux of E 
out of V. With some straightforward calculations, performing an integration exchange justified 
by the uniform convergence already established, one finds: 

jB.ndSo = Jp{P)dVj^^;^dSo+ J p{P)dV J ^f^dSo (19) 

S V 5 ' ' D-V 5 ' ' 

where n is the vector of modulus 1 normal to S. In (19) the second integral does not contain 
singularities, since Pq is on S, while P is external to S. Since a trivial calculation gives 

div(^) = 

when Pq ^ P, a standard application of Gauss Theorem proves that the second integral van- 
ishes. By the same reasoning, in the first integral the boundary S can be modified to a sphere 
centered on P, and totally included in D, without modifying its value. In this conditions, a 
direct calculation shows that 

/ B.ndSo = 47r / p{P)dV 
s V 

and using again Gauss Theorem 
JdivEdF = 47r / p{P)dV 

V V 

Finally, keeping into account that the domain V is arbitrary and that the functions under 
integral are continuous, one obtains 

divE = A-rrp 

which, inserting the definition of E, is equivalent to the equation to be verified (16). The focal 
point in the reported prove is that a vector function like (18) has flux zero from any closed 
surface not including singularities of the integral; from this it follows that the flux from an 
arbitrary surface can be calculated by modifying the surface up to a sphere concentrated on 
the unique singular point, for which the calculation is trivial. Similar methods are applied to the 
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calculation of the residues of an analytical function, providing Cauchy formula, whose aspect 
is similar to (18). This localisation property allows also to understand how the restriction on 
the finiteness of the domain D could be lowered: it is sufficient to separate D into a finite 
part surrounding the point P and a remaining infinite part. For the finite part the reported 
reasoning applies, while the remaining one gives no contribution. Obviously some conditions 
must be added, which authorise the performed interchange of integrals; a sufficient condition 
is the existence of the integrals appearing in (19) as multiple integrals, by Fubini Theorem. 

3.2 Magnetostatic case 

This is the case in which both the density p and the current density j do not depend on time. 
In this case it will be necessary to use condition (5), which becomes 

divj = (20) 

In this case the proposed solutions (14) and (15) assume the form 

Hp.) = -\r-^dv (21) 

D 

v{Po) = -I^^dV (22) 

D 

where D is the finite domain in which p and j do not vanish, while the equations to be verified 
((8), (9) and (10)) reduce to 

divA = (23) 
AA = -^j (24) 
Aif = -Aixp (25) 
(25) has been already proved in the preceding section. The same reasoning, applied to the 
individual components of A and j, proves (24). Therefore all what is required is to prove (23). 
To this aim, the first point is to consider two vectors rg and r and to evaluate the flux of the 
vector from the boundary S of D. Considering firstly the case in which Pq is external 

to D, which avoids any singularity, and using (20), (21), (22) and Gauss Theorem, one obtains 



/ = /div^dy = / ^divjdy + Jj.grad^dy = / -^^dV = divoA(ro) 

S D D D D ' ' 

(in the last passage a derivative under the sign has been taken, justified as in the preceding 
section) . 

Therefore in this case the proof of (23) is obtained, if the flux appearing at the left hand 
side is zero. A sufficient condition for that is that at each point of 5 j is tangent to S, or 
in particular null. To prove this, suppose that exists a point P on 5" at which j is directed 
e.g. towards the extern. Then it should be possible to find an neighborhood of P satisfying 
the same condition (remember that j is assumed continuous), with the consequence that the 
flux of j from such neighborhood, completed to a closed surface, should be strictly positive, in 
contradiction with (20). 

Considering now the case in which Pq is internal to D, the integral appearing in (21) 
becomes singular, and the above passages are not valid. But they can be applied excluding 
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from the integral the contribution of a sphere of ray e centered on ro and with boundary 
Se, finding 

S S Se Se D-De Se 

The first integral is not singular, and can be treated as in the first case; for the second one 
the mean value theorem is applied, obtaining 

where K* is a suitable value comprised between the minimum and the maximum of j.n on 
S, and therefore bounded. Taking the limit for £ — > 0, the desired result divoA(ro) = is 
obtained. 

Finally, in the case of P on the boundary of D, it is suflficient to use the already obtained 
results and the continuity of divA, whose prove is obtained by the method of the preceding 
section. 

3.3 Monocromatic case 

This is the case in which the density p and the current density j have a sinusoidal time de- 
pendence. Using the complex exponential notation, which simplifies some arguments, they are 
given by 



p{x, y, z, t) = pa{x, y, z) exp{-iujt) (26) 
y, z, t) = 3a{x, y, z) exp(-ia;i) (27) 

where Pa{x, y, z) and ia{x, y, z) depend only on the spatial variables, and a; is a fixed parameter. 
In this case the proposed solutions (14) and (15) assume the form 



A(Po) = -lJ ^(-Mp&dV. expi-iut) (28) 
^^p^^ M^,y^W(lMdy. eM-i^t) (29) 

D 

where D is a finite domain outside of which p and j vanish. Posing 
A = Aa exp^—iut) 

if = If a exp(— itui) 

equations (28) and (29) give 



A„(Po) = -i / iA-Mp(mdV (30) 



c 
D 



ip,{p,) = - J e^^^2Mf}p&dv (31) 

D 

while the equations to be verified ((8), (9)and (10)) reduce to 
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divAa - ikcpa = (32) 
A(pa + k'^(pa = -47r99a (33) 

AAa + PAa = -^ia (34) 

where k = w/c.Finally the continuity condition (5) becomes 

divja - iujpa = (35) 
The first step is to verify (32). The detailed steps are not reported, since they are identical 
to those used in the preceding section, with the only difference that the starting point is the 
calculation of the flux of j°'^^pMro-r|) -^^j^if,]^ vanishes due to the tangent condition on ja. 

Coming to the proof of (33), one defines 



Ea = -grad(^a + ikAa (36) 
and, using (32), obtains 



divoEa = -divo /pa(«^|ro - r| - 1) exp(iA;|ro - r\) ^^°_^^.^ dV + k'^tpa 

Using the fact that an integral with first order singularity in [rg — r| can be differentiated 
under the sign and developing some calculations, the following is obtained: 



divoEa = -ikj pa{ik\rQ - r| - l)exp(zA;|ro - v\) ^^^\^^^ dV + divoC (37) 
having defined 



C = /paexp(zfc|ro-r|)^i5^dy 

D 



At this point the flux of C through a surface S contained in V is evaluated by steps similar 
to those reported in the section on electrostatics. The result is 

/ divCdFo = !dS^j Pa{P)^^^ exp(zfc|ro - v\)dV = J+ J Pa{P)dV J ^f-i^ exp(zfc|ro - r\)dSo 

V 5 D ' ' V D-V S 

In the above equation the second integral is regular, and its value is 

/ paiP)dVj^f^eMik\ro-r\)dSo = JdVo J pa(P)^-i^ exp(ifc|ro - r|)dF 



On the contrary, in the first integral the surface is moved up to the already used sphere of 
radius £, and the following is obtained 



/divCdVb = /dyo / p /^-PW-o--l) dy+ / dVo J Pa 7(^1-0-1) + 4nJ padV 

V V D-V V-Ve V V 
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Taking the limit e ^ and using the fact that V is arbitrary, the final result is obtained 

divC = ikj p„22Pg^dy + 47rv.„ 

which inserted into (37) provides the desired result (33). 

The final step is the verification of (34). At a first sight, it could appear sufficient to apply 
the last obtained result to the three components of (34). However, there is a problem, due to 
the fact that to obtain (33) from (31), essential use has been made of (32), which depends on 
the fact that a-priori pa is not arbitrary, but tied to a certain vector by (35). Therefore the 
proof is complete if the fact can be proved that to each component ja' of a vector ja another 
vector j* can be associated, satisfying the equation 



divj* - iuj^, = (38) 

and moreover tangent, as ja, on the boundary of the domain D. 
This is consequence of the following property: 
The equation in j 



divj = /(x,y,2;) (39) 

has solutions, even if j is constrained to be tangent to an assigned surface S. In fact it is 
easy to verify that (39) (but in general not the constraint) is satisfied by the vector 



r = {0,0,1 f{x,y,z)dz) 
and therefore also by 

j = j** + rotP 

where P is an arbitrary vector. Therefore all is reduced to choice P in such a way that on 
rotP.n = — j**.n 

i.e. to prove that the equation in P 

rotP. n = g{x,y, z) (40) 
has solutions on S. 

In fact this is possible even by imposing = Py = 0, and reducing (40) to the form 

dPz dP, . , 



S: 
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which is a first order partial differential linear equation and can be solved through a method 
which transform it into an ordinary differential equation (see any text of Mathematical Analysis, 
e.g. [5]). 



3.4 General case 

This is the case in which the density p and the current density j have arbitrary time dependence, 
but with the restriction to admit a Fourier transform, producing functions vanishing outside 
of a finite domain. Under these conditions, p and j are of the following form 

p{r,t) = -k= J Pa{r,Lo)ex.p{-iu;t)dw (41) 
Kr,t) = ^ Ija{r,uj)exp{-iujt)duj (42) 

where B is a finite domain. Inserting (41) and (42) into (14) and (15) and interchanging 
bounded integrations, the proposed solutions take the form 

A(Po,0 = J exp(-za;t)dn / ■'"('•'^^ exp(iu,d/c) = J A„(r, w) exp(-ia;t)dw (43) 

^(Po,t) = --^JeM-iiot)duJ P'^(r,u>)eMiu>d/c) ^y = 1 J ^^^r, u) eM-i^t)d^ (44) 

B D ^ B 

where and (pa are given by (28) and (29). 

Similarly, taking derivatives under the sign of the regular range bounded integrals (41) 
and (42) and using inversion properties for the Fourier transform, the continuity equation (5) 
provides (32) for A^ and ipa- Using the results of the preceding section, it can be concluded 
that Aa and ipa satisfy (32), (33), (34). Using this fact, one obtains, performing derivatives 
under bounded regular integrals: 



divA + = ^ / (dWAa-'fifa) eM-i^t)di^ = 

AA - ^0 = / (AA„ + ^ A„) eM-i^t)<kv = / j„ exp(-^a;^)da; = j 

D D 

Aip - = f {Aipa + ^ipa) exp{-iLot)duj = J exp{-iLot)du = -Airp 



i.e. (8), (9) and (10) have been proved. 

The conditions of finitely extended Fourier transform on p and j have been imposed since 
they are the simplest way to justify the last passages of interchanging integrations and perform- 
ing derivatives under integration. However, the same justifications are ensured more generally 
by conditions of existence and uniform absolute convergence of multiple integrals of the form 



B D 
B D 

having considered that the factor exp{—iujt), of modulus 1, does not disturb such absolute 
convergence. 
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4 Conclusions 



A proof has been given that the fields calculated applying (6) and (7) to the delayed potentials 
(14) and (15) provide effectively a solution to the Maxwell equations (1), (2), (3), (4). The 
complete details of the proof have been reported using the assumption that the sources p and 
j are regular functions up to the first order derivatives, that they are vanishing outside some 
finite spatial domain, and that their temporal Fourier transform is vanishing outside some 
finite frequency interval. It has also been shown that the last two conditions can be lowered, 
assuming only the multiple summability of certain space/frequency functions constructed from 
the sources. 
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